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Abstract. There has been quite some activity and progress concerning spec- 
tral asymptotics of random operators that are defined on percolation sub- 
graphs of different types of graphs. In this short survey we record some of 
these results and explain the necessary background coming from different 
areas in mathematics: graph theory, group theory, probability theory and 
random operators. 
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1. Preliminaries 

Here we record basic notions, mostly to fix notation. Since this survey is meant 
to be readable by experts from different communities, this will lead to the effect 
that many readers might find parts of the material in this section pretty trivial - 
never mind. 

1.1. Graphs 

A graph is a pair G = (V, E) consisting of a countable set of vertices V together 
with a set E of edges. Since we consider undirected graphs without loops, edges 
can and will be regarded as subsets e = {x,y} C V. In this case we say that e is 
an edge between x and y, respectively adjacent to x and y. Sometimes we write 
x <~ y to indicate that {x,y} <G E. The degree, the number of edges adjacent to x, 
is denoted by 

deg G := deg : V -> N , deg(x) := #{y eV \ x~y}. 

A graph with constant degree equal to k is called a fc-regular graph. 

A path is a finite family 7 := (ei, e2, e„) of consecutive edges, i.e., such 
that efe n ek+\ ^ 0; the set of points visited by 7 is denoted by 7* := e\ U ... U e„. 

To appear in: D. Lenz, F. Sobieczky and W. Woess (Eds.), Boundaries and spectra of random 
walks. Proceedings of the alp workshop Graz - St. Kathrein 2009, Progress in Probability, vol. 
64. Springer, Basel, 2011. 
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This gives a natural notion of clusters or connected components as well as a natural 
distance in the following way. If x is a vertex, then C x , the cluster containing x, 
is the set of all vertices y, for which there is a path 7 joining x and y, i.e., so 
that x, y E 7*. The length of a shortest path joining x and y is called the distance 
dist(x, y). With the convention inf := 00 it is defined on all of V, its restriction 
to any cluster induces a metric. 

A subgraph G' — (V',E') of G is given by a subset V C V and a subset 
E' C £. The subgraph G' = (V, £') induced feu V"' has the edge set = {e e £ | 
e C V'}. 

A one-to-one mapping $ : V — > V is called an automorphism of the graph 
G = (V,E) if {x,y} E E if and only if e E. The set of all auto- 

morphisms Aut(G) is a group, when endowed with the composition of automor- 
phisms as group operation. An action of a group r on G is a group homomorphism 
j : T — > Aut(G), and we write 72; := (i(7))(#) for 7 e T, x E V. An action is called 
/ree, if 72; = x only happens for the neutral element 7 = e of T. A group action 
is called transitive, if the or&zi Tx := {72: | 7 e T} of a; equals V for some (and 
hence every) vertex x E V. Note that in this case G looks the same everywhere. 

Example. A prototypical example is given by the d-dimensional integer lattice 
graph L d with vertex set Z d and edge set given by all unordered pairs of vertices 
with Euclidean distance one. Clearly, the additive group 7L d acts transitively and 
freely on h d by translations. 

For any group action, due to the group structure of T, it is clear that two 
orbits Fx 7^ Ty must be disjoint. If there are only a finite number of different 
orbits under the action of T, the action is called quasi-transitive, in which case 
there are only finitely many different ways in what the graph can look like locally. 
For quasi-transitive actions, there are finite minimal subsets T of V so that 

|J Tx = V. (1.1) 

These are called fundamental domains. 

1.2. The adjacency operator and Laplacians 

The adjacency operator of a given graph G = (V, E) acts on the Hilbert space 
£ 2 {V) of complex- valued, square-summable functions on V and is given by 

A := A G : £ 2 (V) -> £ 2 (V),Af(x) := ]T f(y) for / E £ 2 (V), x E V. 

We will assume throughout that the degree deg is a bounded function on V, and so 
A is a bounded linear operator. The (combinatorial or graph) Laplacian is defined 
as 

A := A G : £ 2 (V) -+ ( 2 {V), A.f(x) := ^[/(x) - f(y)} for / E £ 2 (V), x E V 
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so that A G = D G — A G , where D := Dq denotes the bounded multiplication 
operator with deg. Signs are a notorious issue here: note that (contrary to the 
convention in most of the second author's papers) there is no minus sign in front 
of the triangle. 

For a subgraph G' = (V',E') of a given graph, certain variants of Aqi are 
often considered: The Neumann Laplacian is just Ag, := Ac, meaning that the 
ambient larger graph plays no role at all. The Dirichlet Laplacian Ag, (the nota- 
tion agrees with that of [361 El 13 El]) penalises boundary vertices of G' in G, that 
is vertices with a lower degree in G' than in G: 

Ag, := 2(D G -D G ,)+ Ag, = 2D G - D G , - A G , : l 2 (V) -+ f(V). 

A third variant is called pseudo-Dirichlet Laplacian in [3111 \5T\ ; here we use the 
notation from [5, where it is named adjacency Laplacian: 

A G , := Dq - Dqi + Ag, = D G - A G , : f(V) -► £ 2 {V). 

The motivation and origin for the terminology of the different boundary conditions 
are discussed in - together with some basic properties of these operators. Most 
importantly, they are ordered in the sense of quadratic forms 

< A& < A G , ^ Ag, < 2D G < 2||deg G || 00 Id (1.2) 

on £ 2 (V'). Here, Id stands for the identity operator. We recall that for bounded 
operators on a Hubert space T~L, the partial ordering A ^ B means (ip, (B — 
A)ip) ^ for all ip € %, where the brackets denote the scalar product on H. 
Thus the spectrum of each Laplacian A^, , X £ {N, A,D}, is confined according 
to spec(Ag,) C [0, 2|| deg G ||oo] . The names Dirichlet and Neumann are chosen in 
reminiscence of the different boundary conditions of Laplacians on open subsets of 
Euclidean space. In fact one can easily check that for disjoint subgraphs G\, G2 C 
G, 

A Gi © A G2 < Aq iUG2 < Ag iUG2 < A Gl © A G2 . 

The adjacency Laplacian does not possess such a monotonicity. 

On bipartite graphs, such as the lattice graph h d , the different Laplacians are 
related to each other by a special unitary transformation on £ 2 (V). We recall that 
a graph is bipartite if its vertex set can be decomposed into two disjoint subsets 
V± so that no edge joins two vertices within the same subset. Define a unitary 
involution U = U* = If- 1 on 1 2 {V) by (Uf)(x) := ±f{x) for x £ V±. Clearly, we 
have U*DU = D and U*AU = -A. The latter holds because of 



(A(Uf))(x) = ^(Uf)(y) = £>/(!/) - -(U(Af))(x) 

y^x y^x 
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Figure 1. Two Cayley graphs of Z d . 
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for every x G V±. In particular, for any subgraph G of a fc-regular bipartite graph 
G we get 



(1.3) 



Consequently, spectral properties of the different Laplacians at zero - the small- 
est possible spectral value as allowed by (1.2 1 - can be translated into spectral 
properties (of another Laplacian) at 2k. 

1.3. Amenable groups and their Cayley graphs 

Here we record several basic notions and results that will be used later on; we 
largely follow [5]. 

Let r be a finitely generated group and ScTa symmetric (i.e. S^ 1 C S) 
finite set of generators that does not contain the identity element e of T. The 
Cayley graph G = G(T, S) has T as a vertex set and an edge connecting x, y G T 
provided xy~ x G S. By symmetry of S we get an undirected graph in this fashion, 
and G is | SI -regular. Moreover, it is clear that T acts transitively and freely on G 
by left multiplication. 

Examples. (1) The d-dimensional integer lattice graph L d is the Cayley graph 
of the group Z d (written additively, of course) with the set of generators S = 



{< 



31 



j = l,...,d} with ej the unit vector in direction j. 



(2) Changing the set of generators to S' := S U {±ej ± | 1 ^ j < k ^ d} 
gives additional diagonal edges; see Figure 1 for an illustration in d = 2. 

(3) The Cayley graph of the free group with n G N \ {1} generators gi, g n 
can be formed with S = {g\, g n , , g^ 1 }', it is a 2n-regular rooted infinite 
tree. More generally, a (n + l)-regular rooted infinite tree, k G N \ {1}, is also 
called Bethe lattice M K , honouring Bethe [llj who introduced them as a popular 
model of statistical physics. Every vertex other than the root e in B K possesses 
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one edge leading "towards" the root and n "outgoing" edges, see Figure 2 for an 
illustration for n = 2, respectively K = 3. 

Due to fundamental theorems of Bass [10], Gromov [26] and van den Dries 
and Wilkie [S3], the volume, i.e. the number of elements, of the ball B(n) consisting 
of all those vertices that are at distance at most n from the identity e, 

V(n) := \B(n)\ := #{x e T | dist G(r , s) (x, e) < n}, (1.4) 

has an asymptotic behaviour that obeys one of the following alternatives: 

Theorem 1.1. Let G = G(T, S) be the Cayley graph of a finitely generated group. 
Then exactly one of the following is true: 

(a) G has polynomial growth, i.e., V(n) ~ n d for some d GN. 

(b) G has superpolynomial growth, i.e., for all d G N and b € K there are only 
finitely many n e N so that V(n) ^ bn d . 

The growth behavior, in particular the exponent d, is independent of the chosen 
set S of generators. 

There is another issue of importance to us, amenability. A definition in line 
with our subject matter here goes as follows: 

Definition 1.2. A discrete group T is called amenable, if there is a F0lner sequence, 
i.e., a sequence [F n ) n ^ of finite subsets which exhausts T with the property that 
for every finite F dT: 

\(F-F n )AF n \ 

— -> for n -> oo, 

l-f'nl 

where AAB := (A \ 5) U (B \ A) denotes the symmetric difference of two sets A 
and B. 

There is quite a number of different equivalent characterisations of amenabil- 
ity. The notion goes back to John von Neumann [63] . In its original form he required 
the existence of a mean on £°°(T), i.e., a positive, normed, T-invariant functional. 

Remarks 1.3. (1) The defining property of a F0mer sequence is that the vol- 
ume of the boundary of F n becomes small with respect to the volume of F n itself 
as n — > oo. Boundary as a topological term is of no use here; instead, thinking of 
the associated Cayley graph, F ■ F n can be thought of as a neighborhood around 
F n (at least for F containing the identity) and so \(F ■ F n )AF n \ represents the 
volume of a boundary layer around F n . Thinking of F as the ball B(r) makes this 
picture quite suggestive. 

(2) Discrete groups of subexponential growth are amenable. 

(3) The lamplighter groups (see below) are amenable but not of subexponential 
growth. Consequently, growth does not determine amenability. 

(4) The standard example of a nonamenable group is the free group on two 
generators. 
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Figure 2. Bethe lattice B3, the Cay ley graph of the free group 
with n — 2 generators a, b. 

Let us end this subsection with the example we already referred to above: 
Example. Fix m£N,m)2. The wreath product Z m ; Z is the set 
Z m ; Z := {(if, x) I </? : Z — > Z m , supp </? finite , a; € Z}, 

(Vl,3fl) * (f2,X 2 ) := {<Pl +tp2(-- Xi),Xi + X 2 ) 

and is called the lamplighter group. It is amenable, see [7]. 
2. Spectral asymptotics of percolation graphs 

This section contains the heart of the matter of the present survey. After intro- 
ducing percolation, we begin discussing the relevant properties of the random op- 
erators associated with percolation subgraphs. The central notion is the integrated 
density of states, a real- valued function. We then explain a number of results on the 
asymptotic behaviour of this function and how methods from analysis, geometry 
of groups, graph theory and probability are used to derive these results. 

2.1. Percolation 

Percolation is a probabilistic concept with a wide range of applications, usually 
related to some notion of conductivity or connectedness. Its importance in (statis- 
tical) physics lies in the fact that, despite its simplicity, percolation yet exposes a 
phase transition. The mathematical origin of percolation can be traced back to a 
question of Broadbent that was taken up in two fundamental papers by Broadbent 
and Hammersley in 1957 15, 28J. Percolation theory still has an impressive list of 
easy-to-state open problems to offer, some with well established numerical data 
and conjectures based on physical reasoning. We refer to [UJ EH] f° r standard 
references concerning the mathematics, as well as Kesten's recent article in the 
Notices of the AMS [32]. 

Mathematically speaking, and presented in accordance with our subject mat- 
ter here, percolation theory deals with random subgraphs of a given graph G = 
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Figure 3. Part of a realisation G w for bond percolation on L 2 
for p = \. 

(V, E) that is assumed to be infinite and connected. A good and important example 
is the d-dimensional lattice graph h d , the particular case d = 1 being very special, 
however. There are two different but related random procedures to delete edges 
and vertices from G, called site percolation and bond percolation. In both cases, 
everything will depend upon one parameter p € [0, 1] that gives the probability of 
keeping vertices or edges, respectively. 

Let us start to describe site percolation. We consider the infinite product 

n := n s[tc := {0, l} y , P p := (g) (p • S 1 + (1 - p) ■ 5 ) , 

x£V 

as probability space with elementary events lu := (uj x ) xe v, w x G {0,1}, and a 
product Bernoulli measure P p that formalizes the following random procedure. 
Independently for all vertices (also called sites in this context) of V, we delete the 
vertex x from the graph with probability 1 — p, along with all edges adjacent to x. 
This corresponds to the event cu x = 0, and we call the site x closed. On the other 
hand, we keep the vertex x and its adjacent edges in the graph with probability 
p. This corresponds to the event uj x = 1, in which case we speak of an open 
site. Every possible realisation or configuration is given by exactly one element 
uj = (uj x ) X £v € 51, and the measure ¥ p above governs the statistics according to 
the rule we just mentioned. Note that we omit the superscript in the notation of 
the product measure. The graph we just described is illustrated in Figure 3 and 
formally defined by G u = {V w , E w ), where 

V u := {x e V | lj x = 1}, E UJ :={eeE\e<Z V u }, 

i.e. the subgraph of G induced by V u . Note that for p = the graph G w is empty 
with probability 1 and for p — 1 we get G^ — G with probability 1. 
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The second variant, bond percolation, works quite similarly: 

Q := n hoad := {0, 1} E , F p := ■ S x + (1 - p) ■ 5 ), 

leading to the subgraph G w = (V W ,E U ) with 

V u := V, := {e e E \ u e = 1}. 

It amounts to deleting edges (also called bonds in this context) with probability 
1 — p, independently of each other. The choice V u — V is merely a convention. 
Other authors keep only those vertices that are adjacent to some edge. 

In both site and bond percolation, the issue is the connectedness of the so- 
obtained random subgraphs. Note that the realisations G u themselves do not de- 
pend upon p, while assertions concerning the probability of certain events or the 
stochastic expectation of random variables constructed from the subgraphs surely 
do. A typical question is whether the cluster C x that contains vertex x € V is 
finite in the subgraph G u for P p -almost all uj £ 51 or whether it is infinite with 
non-zero probability. In the latter case one says that percolation occurs. 

Let us assume from now on that G is quasi-transitive, so that the above 
question will have an answer that is independent of x. The percolation threshold 
or critical probability is then defined as 

p H :=sup{pe [0,1] | F p [\C x \ = oo] = 0}. 



It is independent of x since, globally, G looks the same everywhere, cf. (1.1), and 
P p is a product measure consisting of identical factors. A related critical value is 
given by 

p T :=sup{pe [0,1] | Ep[|C x |] <oo}, 
and it is clear that pt ^ ph- Here, E p stands for the expectation on the probability 
space (fi,P p ). The equality of these two critical values is often dubbed sharpness 
of the phase transition, and we write p c := pn = Pt in this case for the critical 
probability. Clearly, sharpness of the transition is a desirable property, as both 
Ph and pt represent two equally reasonable ways to distinguish a phase with P p - 
almost surely only finite clusters, the subcritical or non-percolating phase, from a 
phase where there exists an infinite cluster with probability one, the supercritical 
or percolating phase. Apart from that, sharpness of the phase transition has been 
used as an important ingredient in the proof of Kesten's classical result that p c — 
| for bond percolation on the 2-dimensional integer lattice L 2 . Together with 
estimates known for p < px, it gives that the expectation of the cluster size decays 
exponentially, i.e., 

V P {\C x \=n} ^e- a - n , neN, 
with some constant a p > for all p < p c . This fact is also heavily used in some 
proofs of Lifshits tails for percolation subgraphs, see below. Fundamental papers 
that settle sharpness of the phase transition for lattices and certain quasi-transitive 
percolation models arc 2, 47, 48 . Recent results valid for all quasi-transitive graphs 
can be found in [6] together with a discussion of the generality of earlier literature. 
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Theorem 2.1. ([B], Theorem 2, Theorem 3) For every quasi-transitive graph 

PT =PH =-Pc, 

and for every p < p c there exists a constant a p > so that 

P P {\C X \ ^ n} < e~ a " n for allx£V,ne N. 

It is expected that sharpness of the phase transition also holds for perco- 
lation on more general well-behaved graphs even without quasi-transitivity. The 
celebrated Penrose tiling gives rise to such a graph without quasi-transitivity but 
some form of aperiodic order. A result analogous to Thm. |2.1| was proven for the 
Penrose tiling in [3D]. The general case of graphs with aperiodic order has not yet 
been settled. We refer to [3S] for partial results in this direction. 

2.2. The integrated density of states 

The study of the random family (A(j ij) ) UJe n of Laplacians on percolation graphs 
was proposed by de Gennes [HO [2D] and often runs under the header quantum 
percolation in physics. In this paper we focus on the integrated density of states 
(IDS), also called spectral distribution function, of this family of operators. 

In general, the IDS is the distribution function of a (not necessarily finite) 
measure on R that is meant to describe the density of spectral values of a given self- 
adjoint operator. In the cases of interest to us here, the underlying Hilbert space is 
£ 2 (V), with V being the countable vertex set of some graph. In this situation the 
IDS is even the distribution function of a probability measure on R, as we shall see. 
Before giving the rigorous definition that applies in this setting, let us first start 
with a discussion at a heuristic level. For elliptic operators acting on functions 
on some infinite configuration space V with a periodic geometric structure, one 
typically does not have eigenvalues, but rather continuous spectrum. However, 
the restrictions of these operators to compact subsets K of configuration space 
V (more precisely to P(K), actually) come with discrete spectrum. Therefore, 
one can count eigenvalues, including their multiplicities. The idea of the IDS is 
to calculate the number of eigenvalues per unit volume for an increasing sequence 
K n of compact subsets and take the limit. For this procedure to make sense, the 
operator has to be homogenous, at least on a statistical level. Two situations 
are typical: Firstly, a periodic operator, quite often the Laplacian of a periodic 
geometry. And, secondly, an ergodic (statistically homogenous) random family of 
operators, in which case the above mentioned limit will exist with probability one. 

Let H be a self-adjoint operator in £ 2 (V). An intuitive ansatz for the definition 
of the IDS might be N : R -> [0, 1], 



\F n \ n^oo \F n \ 

. (2 ' 1} 

where (F n ) nl zfi is an appropriate sequence of finite sets exhausting V. Before we 



go on, let us add some remarks on our notation in (2.1). In general, we write 1a 



for the indicator function of some set A. Above, 1^ is to be interpreted as the 
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multiplication operator corresponding to the indicator function 1 p n ■ In view of 
the functional calculus for self-adjoint operators we write 1b(H) for the spectral 
projection of H associated to some Borel set B CI. Finally, tr stands for the 
trace on £ 2 (V) and 5 X G £ 2 (V) for the canonical basis vector that is one at vertex 
x and zero everywhere else. 

As was already mentioned, a certain homogeneity property is necessary in 



order for the limit in (2.1| to exist. A careful choice of the exhausting sequence 
is necessary, too. For amenable groups tempered F0lner sequences will do the 
job, as is ensured by a general ergodic theorem of Lindenstrauss [43]. We refer to 
|39l 1491 150] for more details in the present context and sum up the main points in 
the following definition and the subsequent results. 

Definition 2.2. Let G be a graph and let T be an infinite group that acts quasi- 
transitively on G. We fix a fundamental domain T . For EeMwe define 

N peI (E) := ± tr [l^-^E] (A g )] (2.2) 

to be the IDS of the full graph. Secondly, the expression 

N X (E) := N { ^(E) := ± E p {tr [l^]-oo, E]( A G J]} (2-3) 

is the IDS of the Laplacians on random percolation subgraphs, where X £ {N, A, D} 
stands for one of the possible boundary conditions discussed in Subsection |1.2| 



Remarks 2.3. (1) We could have chosen a more general probability measure 
than P p , as long as it is invariant under T. 

(2) Usually, we will omit the superscript p and write simply Nx for the quan- 



tity in (2.3) 



(3) Note that iV por = for any X £ {N, A, D}. 

(4) Note also that Nx is not defined in terms of a single operator , but 
rather using the whole family (Ag )uien'i see also the subsequent result for a 
clarification. 

The next theorem establishes the connection between the heuristic picture 



displayed in (2.1 1 and the preceding definition. The point here is the generality 
of the group involved. In the more conventional setting of random operators on 
Euclidean space M. d (with the group action of Z d ), the equation is the celebrated 
Pastur-Shubin trace formula. 

Theorem 2.4. ([39 , Theorem 2.4) Let G be a graph and let T be an infinite group 
that acts quasi-transitively on G. Then there is a sequence (F n ) ne fa of finite subsets 
of V so that 

N X (E) = lim p^tr [l,^,^ (1 F „ A* 1 F J] , (2.4) 
uniformly in E el for P p -a.e. well. 
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Remarks 2.5. (1) We refer to [H |M1 13H1 EH HI] for further predecessors of 
the latter theorem. 



(2) The inequalities in (1.2 1 imply 

N D sC N A N N . 

(3) A comprehensive theory of the IDS in the (more conventional) set-up of 
random Schrodinger operators can be found in the monographs [17l [531 I2H] > see 
also the surveys [3H [351 IE2] an d the references therein. 

Interestingly, the IDS links quite a number of different areas in mathematics: 
We started with an elementary operator theoretic point of view. If we rephrase 
the basic existence problem in the way that we regard the counting of eigenval- 
ues as evaluating the trace of the corresponding eigenprojection, we arrive at the 
question, whether appropriate traces exist on certain operator algebras. Typically, 
the operators we have in mind are intimately linked to some geometry, so that 
quantities derived from the IDS play an important role in geometric analysis. An 
important example is the Novikov-Shubin invariant of order zero, which equals the 
van Hove exponent in the mathematical physics language and will be discussed in 
our setting further below; see [53] [57] and the Oberwolfach report [21 . Another 
wellknown principle provides a link to stochastic processes and random walks: The 
Laplace transform of Njy is the return probability of a continuous time random 
walk on the graph; details geared towards the applications we have in mind can 
be found in [5Tj . 

The original motivation and the name IDS come from physics. The Laplacians 
we consider show up as energy operators for a quantum-mechanical particle which 
undergoes a free motion on the vertices of the graph. If v, v' G V are connected 
by an edge, the particle can "hop" directly from v to v' or vice versa. In this 
way, the spectrum of the Laplacian appears as the set of possible energy values 
the particle may attain, hence the name IDS for the quantities in Def. |2.2| In the 
percolation case, the motion is interpreted to be a quantum mechanical motion of 



a particle in a random environment. Thm. 2.4 is interpreted as the self- averaging 
of the IDS for a family of random ergodic operators: for P-a.e. realisation u> of the 
environment, the normalised finite- volume eigenvalue counting function converges 
to a non-random quantity. In particular, if one had taken an expectation on the 



r.h.s. of (2.4), one would have ended up with the very same expression in the 
macroscopic limit. 

The IDS is one of the simplest, but nonetheless physically important spectral 
characteristics of the operators we consider. It encodes all thermostatic properties 
of a corresponding gas of non-interacting particles. As an example we mention 
a systems of electrons in a solid, where this is a reasonable approximation in 
many situations. Besides, the IDS enters transport coefficients such as the electric 
conductivity and determines the ionisation properties of atoms and molecules. For 
this reason, the IDS (more precisely, its derivative with respect to E, the density 
of states) is a widely studied quantity in physics. 
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2.3. The integer lattice 

In this subsection we are concerned with the asymptotics at spectral edges of the 
IDS of the family of Laplacians (Ag ) w ef2 orL bond-percolation subgraphs of the 
d-dimensional integer lattice graph h d (or bond percolation on Z d , for short). 

The spectral edges of these Laplacians turn out to be and 4d. In fact, 
standard arguments which are based on ergodicity w.r.t. Z d -translations, 
yield that even the whole spectrum equals almost surely the one of the Laplacian 
A^d on the full lattice 

spec(Ag ) = [Q, 4d] for P p -almost every ui G f2, 

any p s]0, 1] and X E {N, A, D}. Thus, the left-most and right-most inequality in 



(1.2 1 are sharp in this case. Since the lattice L d is bipartite, it follows from (1.3 1 
with k — 2d that the different Laplacians are related to each other by a unitary 
involution, which implies the symmetries 

N A (E) = 1- lim N A (s), 

£t (2 5) 

N D(N) (E) = 1- lim N N(D) (s) 

e^id—E 

for their integrated densities of states for all E € [0,4d]. The limits on the right- 



hand sides of (2.5) ensure that the discontinuity points oiNx are approached from 
the correct side. 

As before we write p c = p c (d) for the unique critical probability of the bond- 
percolation transition in 7L d . We recall from [53] that p c = 1 for d = 1, otherwise 
p c G]0, 1[. Let us first think about what to expect. At least for small p, the ran- 



dom graph G u is decomposed into relatively small pieces, due to Theorem 2.1 
above. This means that there cannot be many small eigenvalues as the size of 
the components limits the existence of low lying eigenvalues. Consequently, the 
eigenvalue-counting function for small E must be small. It turns out that the IDS 
vanishes even exponentially fast. This striking behaviour is called Lif shits tail, 
to honour Lifshits' fundamental contributions to solid state physics of disordered 
systems [40] [4TJ [42] . In fact, Lifshits tails continue to show up in the percolating 
phase for the adjacency and the Dirichlet Laplacian at the lower spectral edge. 
This follows from a large- deviation principle. 

Theorem 2.6. ([51], Theorem 2.5) Assume d E N andp e]0, 1[. Then the integrated 
density of states Nx of the Laplacians (A^ ) w ga on bond-percolation graphs in 
Z d exhibits a Lifshits tail at the lower spectral edge 

is «a._^ f0TXe{AD} ( , 6) 

and at the upper spectral edge 

hm ^- Nx ™ =-$ for Xe{N,A } . (2.7) 
E^id ln(4d — E) 2 v ' 
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Actually, slightly stronger statements without logarithms are proven in [EH 



see the next lemma. Together with the symmetries (2.5), these bounds will imply 
the above theorem. 

Lemma 2.7. ([SI], Lemma 3.1) For every d £ N and every p £}0, 1[ there exist 
constants Sd, ol u , ai €]0,oo[ such that 

exp{- a/ £T d/2 } sC N D (E) < N A (E) < exp{-a n £T d / 2 } (2.8) 

holds for all E €]0,£_d[. 

Remarks 2.8. (1) In the non-percolating phase, p £]0,p c [, the content of The- 
orem 2.6 has already been known from [36] . where it is proved by a different 
method. The method of [36], however, does not seem to extend to the critical 
point or the percolating phase, p £\p c , 1[. 



(2) The Lifshits asymptotics of Theorem 2.6 are determined by those parts 
of the percolation graphs which contain large, fully-connected cubes. This also 
explains why the spatial dimension enters the Lifshits exponent d/2. 



(3) We expect that (2.6) can be refined in the adjacency case X = A as to 
obtain the constant 

»Bf#=-. CM- <*> 

An analogous statement is known from Thm. 1.3 in [12] for the case of site- 
percolation graphs. Moreover, it is demonstrated in [4] that the bond- and the 
site-percolation cases have similar large-deviation properties. 



The second main result of this subsection complements Theorem |2.6| in the 
non-percolating phase. 

Theorem 2.9. ( 36J, Theorem 1.14) Assume d £ N and p £]Q,p c [. Then the inte- 
grated density of states of the Neumann Laplacians (Ag )^gn on bond-percolation 
graphs in Z d exhibits a Lifshits tail with exponent 1/2 at the lower spectral edge 

lnm^)-^(0)]| = _l 
Eia InE 2 ' v ' 

while that of the Dirichlet Laplacians (A^ )uien exhibits one at the upper spectral 
edge 

\n\\n[Nn(4d)-N D (E)}\ _ 1 
E™d ln(4d-E) ~ 2' 1 ' 

where N D (Ad) := lim Bf4(j N D (E) = 1-N N (0). 

Remarks 2.10. (1) This theorem also follows from sandwich bounds analogous 
to those in Lemma l2~7l We do not state them here but refer to Lemmas 2.7 and 2.9 
in [515] for details. Using interlacing techniques, [SS] establishes a better control on 
the constants in these bounds. For example, it was found that for all sufficiently 
small energies E 

N N (E) - N N (0) < AEexp{-a + £T 1/2 } (2.12) 
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with a+ := 4/[3v3~Xp]j where \v stands for the expected number of vertices in 
the cluster containing the origin and where the constant A > can also be made 
explicit. 

(2) The constant A^v(O) appearing in Theorem |2.9| is given by 



A^at = im ^ r- = p{p) + (1 - p 2.13 

A|Z d |A| 

and equals the mean number density p(p) of clusters with at least two and at most 
finitely many vertices, see e.g. Chap. 4 in [25], plus the number density of isolated 
vertices. This follows from the fact that the operator l[o,oo[( — Ag^ A ) is nothing 
but the projector onto the null space of the restriction A of Ag^ to l 2 (K). The 
dimensionality of this null space equals the number of finite clusters and isolated 
vertices of in A, see Remark 1.5(iii) in [36] . 

(3) The Lifshits tail for Nn at the lower spectral edge - and hence the one 
for Nd at the upper spectral edge - is determined by the linear clusters of bond- 
percolation graphs. This explains why the associated Lifshits exponent —1/2 is 
not affected by the spatial dimension d. Technically, this relies on a Cheeger in- 
equality [18] for the second-lowest Neumann eigenvalue of a connected graph, see 
also Prop. 2.2 in [3"6l . 

The third main result of this subsection is the counterpart of Theorem [279] in 
the percolating phase. 

Theorem 2.11. ([5T]. Theorem 2.7) Assume d E N \ {1} and p e]p c ,l[. Then 
the integrated density of states of the Neumann Laplacians (A.q ) u gs) on bond- 
percolation graphs in 7L d exhibits a van Hove asymptotic at the lower spectral edge 

\n[N N (E) — jVjy(O)] = d 
Eio ]nE 2 ' y ' 

while that of the Dirichlet Laplacian Ag^ exhibits one at the upper spectral edge 

HN D (ld)-N D (E)] d 
E$d ln(4d - E) ~ 2' ( • 5J 



Similar to the two theorems above, Theorem 2.11 also follows from upper 



and lower bounds and the symmetries (2.5). 



Lemma 2.12. ([ST], Lemma 4.1) Assume d g N \ {1} and p £}p c , 1[. Then there 
exist constants En, C u , Ci €]0, oo[ such that 

dE d ' 2 < N N {E) - N N {G) < C u E d ' 2 (2.16) 

holds for all E g]0, £jv[- 



Remarks 2.13. (1) Lemma 2.12 relies mainly on recent random- walk estimates 
[46], El [29] for the long-time decay of the heat kernel of on the infinite cluster. 
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(2) There is also an additional Lifshits-tail behaviour with exponent 1/2 due 
to finite clusters as in Theorem |2.9[ but it is hidden under the dominating van Hove 
asymptotic of Theorem |2.11| Loosely speaking, Theorem |2.11| is true because the 
percolating cluster looks like the full regular lattice on very large length scales 
(bigger than the correlation length) for p > p c . On smaller scales its structure is 
more like that of a jagged fractal. The Neumann Laplacian does not care about 
these small-scale holes, however. All that is needed for the van Hove asymptotic 
to be true is the existence of a suitable ci-dimensional, infinite grid. The adjacency 
and Dirichlet Laplacians though do care about those small-scale holes, as we infer 
from Theorem 12.61 

(3) In the physics literature the terminology van Hove "singularity" is also 
used for this kind of asymptotic. This refers to the fact that for odd dimensions d 
derivatives seize to exist for high enough order. 

The above three theorems cover all cases for p and X except the behaviour at 
the critical point p — p c of iVjy at the lower spectral edge, respectively that of Np 
at the upper spectral edge. In dimension d — 2 upper and lower power-law bounds 
have been obtained in |57j . However, the exponents differ so that the asymptotics 



is still an open problem; see also Remark 2.17(3) below for further properties at 
criticality. 

2.4. The regular infinite tree (Bethe lattice) 

In this subsection we report results from [55 on the asymptotics at spectral edges 
for the IDS of the family of Laplacians ( Ag ) we n on bond-percolation subgraphs of 
the (k + l)-regular rooted infinite tree, a.k.a. Bcthe lattice B K , where k£N \ {1}. 
Percolation on regular trees is well studied, see e.g. [54 , and it turns out that 
the bond-percolation transition occurs sharply at the unique critical probability 



p c = k . Here, sharpness of the phase transition is implied by, e.g., Theorem 2.1 
but it can also be verified by explicit computations. In contrast to percolation on 
the hypercubic lattice h d , where the infinite cluster of the percolating phase is 
unique, there exist infinitely many percolating clusters simultaneously for p > p c 
on B K . 

The results on spectral asymptotics of the IDS are analogous in spirit to 
the ones of the previous subsection, but restricted to the non-percolating phase. 
However, as the Bethe lattice B K exhibits an exponential volume growth of the 
ball B(n) of radius n about its root 

n , 

V(n) = \B(n)\ = 1 + (« + 1) = 1 + («" - 1) !LL T, 

i K — l 

cf. Figure [2j there will be natural differences. 

The next lemma determines the spectral edges of the operators under consid- 
eration. As a consequence of the exponential growth of the graph, and in contrast 
to the preceding subsection, the spectrum of the Laplacian on the Bethe lattice 
does not start at zero, neither does it extend up to twice the degree 2(k + 1). 
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Lemma 2.14. Let k 6 N \ {1} and let Ab k be the Laplacian on the (full) Bethe 
lattice B K . Then 

spec(A B J = [E-,E+], where E% := (^± l) 2 . 

Moreover, for V-almost every realisation G u of bond-percolation subgraphs of B K 
we have 

spec(AgjC[0,£+], spec(A^) = [E~,E+], spec(AgJ C [E~,2(k+1)]. 

Remarks 2.15. (1) We believe that equality (and not only "C") holds for the 
statements involving the Neumann and the Dirichlet Laplacians, too. 

(2) Since the Bethe lattice is bipartite the above lemma reflects the symmetries 

(Of. 

(3) Almost-sure constancy of the spectra (i.e. independence of u>) is again a 
consequence of ergodicity of the operators, see e.g. [I] for a definition of the crgodic 
group action. 

The ergodic group action on the Bethe lattice, which was referred to in the 
last remark above, is even transitive so that the IDS iVjj of the family (Aq ) w en 
can be defined as in Definition |2.2| with the fundamental cell F consisting of just 
the root. Clearly, Nx will then obey the symmetry relations 

N A {E) = l- lim N A {e), 

et2{K+l)-E 

N D{N) (E) = 1- lim N N{D) (e) (2 ' 1?) 

for all E £ [0,2(k + 1)]. 

Our first result concerns the asymptotic of Njy at the lower edge, resp. of iVo 
at the upper edge. Since these two spectral edges are unaffected by the exponential 
volume growth, it comes as no surprise that we find the same type of Lifshits tail 
as in the Z d -case. 

Theorem 2.16. ([SS]) Assume K € N \ {1} and p e]Q,p c [. Then the integrated 
density of states of the Neumann Laplacians (Ag ) w gn on bond-percolation graphs 
in B K exhibits a Lifshits tail with exponent 1/2 at the lower spectral 



ln\HN N (E)-N Nm \ = _l 
eiq InE 2 ' v ' 

while that of the Dirichlet Laplacian Ag exhibits one at the upper spectral edge 

ln|ln[^(2(. + l))-^)]|_l 
st2(«+i) ln(2( K + 1)-E) 2 

where N d {2{k + 1)) := lim BT2(K+1) N D {E) = 1 - N N (0). 

Remarks 2.17. (1) These asymptotics are again determined by the linear clus- 
ters of bond-percolation graphs, cf. Remark 2.10 (3) The interpretation of the ref- 
erence value Nn (0) in terms of the cluster plus isolated vertex density is analogous 



to Remark 2.10(2) 
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(2) In contrast to this Lifshits-tail behaviour in the subcritical phase, one 
expects Nn(E) — Nn(0) to obey a power-law for small E at the critical point p c , 
caused by the finite critical clusters. This is not yet fully confirmed, but upper and 
lower algebraic bounds (with different exponents) follow from the random-walk 
estimates in [57] . 

(3) It should be noted that the power-law behaviour at p c mentioned in the 
previous remark is not the one referred to by the famous Alexander- Orbach con- 
jecture [3]. The latter concerns the £ 4 / 3 -behaviour as E -> of N N (E) on the 
incipient infinite percolation cluster. For the case of the Bethe lattice this asymp- 
totic was proven in [9]. (Here no subtraction of Nn(0) is necessary. Instead, one 
kind of conditions on the event that the origin belongs to an infinite cluster, see e.g. 
P3| for details of the definition.) The Alexander- Orbach conjecture says that the 
_B 4 / 3 -asymptotic should also hold for percolation in Z d for every d ^ 2. Extensive 
numerical simulations indicate that this is not true in d = 2 [53] . We refer to [IB] 
for a comprehensive discussion and further references from a Physics perspective. 

In order to reveal the characteristics of the Bethe lattice we now turn to the 
spectral edges E^. 

Theorem 2.18. ([SS]) Assume K € N \ {1} and p E]0,p c [. Then the integrated 
density of states o/(Ag^) ue n on bond-percolation graphs inM K exhibits a double- 
exponential tail with exponent 1/2 at the lower spectral edge 

In [In I In iV x (£011 1 

lim — l - — 1 V ;lJ = — for X = A,D 2.20 

EiE- HE-Ek) 2 v ' 

and one at the upper spectral edge 

In [In I In (1 - N X (E))\] 1 
lim — i — LV s { >,U = -- forX = N.A. 2.21 



Remarks 2.19. (1) The extremely fast decaying asymptotic of (2.20) - and 
similarly that of (2.21 ) - is determined by the lowest eigenvalues E ~ E~ + R~ 2 



of those clusters in the percolation graph which are large fully connected balls 
of radius R. Their volume is exponentially large in the radius, V(R) ~ e R ~ 
c( E ~ E * ) , and their probabilistic occurrence is exponentially small in the vol- 
ume. 



(2) One would expect Theorem 2.18 to be valid beyond the non-percolating 



phase. However, the region p J? p c is still unexplored. 



(3) A double-exponential tail as in ( 2.20 ) will also be found in Theorem 2.24 (3) 



below. This concerns the lower spectral edge of the IDS for percolation on the Cay- 
ley graph of the lamplighter group, which is amenable. These double-exponential 
tails in two concrete situations should also be compared to the less precise last 
statement of Theorem |2.21| below, which, however, holds for superpolynomially 
growing Cayley graphs of arbitrary, finitely generated, infinite, amenable groups. 
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2.5. Equality and non-equality of Lifshits and van Hove exponents on amenable 
Cayley graphs 

... is almost the title of a paper by Antunovic and Veselic [7]. Here we record their 



main results. In our definition of the IDS in Subsection 2.2 above, two entirely 
different cases were treated. Let us first consider the deterministic case of the 
Laplacian on the full graph, denoted by N pcr . In our case of a quasi-transitive 
graph the geometry looks pretty regular; just like in the case of a lattice, the local 
geometry has the same local structure everywhere. Specializing to Cayley graphs 
this allows one to relate the asymptotic of -ZV por near to the volume growth V(n) 



defined in (1.4). The latter is the same for the different Cayley graphs of the same 



group, see Theorem |1. 1| above. 

Theorem 2.20. Let T be an infinite, finitely generated, amenable group, G = 
G(r, S) a Cayley graph of T and N per the associated IDS. If G has polynomial 
growth of order d, then 

IniVperCE) = d 

hxE 2' 



lim 



(2.22) 



// G has superpolynomial growth, then 



lim 



In N peT (E) 
InE 



Proofs can be found in [60l [44] . Note that the limit appearing in (2.22) is 



exactly the zero order Novikov-Shubin invariant, where zero order refers to the 
fact that we deal with the Laplacian on 0- forms, i.e., functions. 

Next we turn to the asymptotic of the IDS Nx of the corresponding perco- 
lation subgraphs. Again, Lifshits tails are found. 

Theorem 2.21. ([7], Theorem 6) Let G = G(T,S) be the Cayley graph of an in- 
finite, finitely generated, amenable group. Let Nx be the IDS for the Laplacians 
)uen of percolation subgraphs of G with boundary condition X G {A, D} in 



the subcritical phase, i.e., for p < p c . 
all E > small enough 

N D (E) < N A {E) ^ exp 



Then there is a constant a p > so that for 



2 



2 V / 2|S*| 



E' 



1 



where V(t) :— V([t\), the volume V(n) is given by ( |1.4[ ) and \ t\ denotes the integer 
part of t G K. If G has polynomial growth of order d, then there are constants 
aj, cvq > so that for E > small enough 



exp 



^N D (E)^N A {E)^exp 



-a 



D 



E- 



If G has superpolynomial growth, then 
In | \nN D (E)\ 



lim 



]nE\ 



lim 

ElO 



]n\ln N A (E)\ 
I In .El 



(2.23) 
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Theorem 2.18 and Theorem 2.24 (3) provide much more detailed information 
as compared to (2.23), but only in two specific situations: the non-amenable free 



group with n 2 generators and the amenable lamplighter group. 

The equality that is mentioned in the title of this subsection is now an easy 
consequence. 

Corollary 2.22. In the situation of the preceding theorem the van Hove exponent 
and Lifshits exponents for X € {A, D} coincide, i.e., 



lim 

ElO 



\\nE\ 



lim 

BJ.0 



lnllnAUQE)! 
I In SI 



lim^H 
Eio m E 



.(E) 



Note that the asymptotic proved for Njj and Na in the case of polynomially 
growing Cayley graphs is actually more precise than the double-log-limit that ap- 
pears in the preceding corollary. For Cayley graphs with superpolynomial growth, 
a lower estimate is missing. However, for the lamplighter groups a more precise 



statement can be proven, see Theorem 2.24 below 



The results of the previous section for the lattice case indicate that one should 
expect a different behaviour for the IDS Njy of the Neumann Laplacian at the 
lower spectral edge: it should be dominated by the linear clusters for p < p c . This 
is indeed true. 

Theorem 2.23. ([7], Theorem 14) In the situation of the previous theorem there 
exist constants a^,ajj > so that for all E > small enough 



exp 



-a N E 2 



^ N N {E) - N N {0) <exp 



-a+E-z 



The dimension d is replaced by 1 in these estimates, since linear clusters are 
effectively one-dimensional and independent of the volume growth of G. This latter 
result remains true for quasi-transitive graphs with bounded vertex degree. 

As already announced, here are the more detailed estimates for the lamp- 
lighter group. 

Theorem 2.24. ([7j, Theorems 11 and 12) Let G be a Cayley graph of the lamp- 
lighter group Z m I 7L. 

(1) There are constants a^,ajj~ > so that for all E > small enough 



N peT (E) sC a^exp 



(2) For every r > | there are constants a l r ,a 2 _ r > so that for all E > 
small enough 

N pcl .(E) ^ a- r exp [-a^ r E~i] . 

(3) For every p < p c there are constants &i, &2> Ci, C2 > so that for all E > 
small enough 



exp 



-c\e 



C 2 E 2 



s? N D {E) < N A {E) ^ exp 



-he 



b 2 E~2 
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2.6. Outlook: some further models 

To conclude, we briefly mention two other percolation graph models for which 
the Neumann Laplacian exhibits a Lifshits-tail behaviour with Lifshits exponent 
| at the lower spectral edge E — in the non-percolating phase. As in the cases 
we discussed above, see Theorem |2.9| for the integer lattice, Theorem 1 2 . 1 6| for the 
Bethe lattice and Theorem |2.23| for amenable Cayley graphs, these Lifshits tails 
will also be caused by the dominant contribution of linear clusters. For this reason 
they occur quite universally, as long as the cluster-size distribution of percolation 
follows an exponential decay - no matter how complicated the "full" graph G may 
look like. This structure will not be seen by the linear clusters of percolation! 

The first class of models [50] consists of graphs G which are embedded 
into M. d (or, more generally, into a suitable locally compact, complete metric space) 
with some form of aperiodic order. The celebrated Penrose tiling in IR 2 constitutes 
a prime example. But one can consider rather general graphs whose vertices form 
a uniformly discrete set in M. d and whose edges do not extend over arbitrarily 
long distances. Amazingly, the main point that needs to be dealt with to establish 
Lifshits tails for such models concerns the definition of the IDS. In contrast to 



the definition in (2.3), one cannot expect to benefit from a quasi-transitive group 
action on G with a finite fundamental cell in this aperiodic situation. The way out 
is to consider the hull of the graph G, that is the set of all Retranslates of G, closed 
in a suitable topology which renders the hull a compact dynamical system. As such 
it carries at least one K^-ergodic probability measure /i, and the expectation in 



(2.3 1 will be replaced by a two-stage expectation: one with respect to /i over all 



graphs G' in the hull of G, and inside of it, for each graph G', the expectation 
over all realisations of percolation subgraphs of G'. The interested reader is 



referred to (37] [50] for more details. 

The second model, Erdos-Renyi random graphs \23\ I13j. has a combinatorial 
background. There we consider bond percolation on the complete graph K n over 
n vertices with bond probability p :~ c/n. The n-independent parameter c > 
corresponds to twice the expected number density of bonds, if n is large. This 
is sometimes referred to as the (very) sparse case. For c €]0, 1[, the fraction of 
vertices belonging to tree clusters tends to 1 as n — > oo, and the limiting cluster- 
size distribution decays exponentially. In this model the IDS is defined by 

N N (E);= lim E<£> [{8i, l].^] (A& )5 1 )] , 
and it exhibits a Lifshits tail at the lower spectral edge E = with exponent 1/2 
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